Renewable  and  Sustainable  Energy  Reviews  13  (2009)  2012-2021 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Renewable  and  Sustainable  Energy  Reviews 

journal  homepage:  www.elsevier.com/locate/rser 


Reliability  and  cost  analyses  of  electricity  collection  systems  of  a  marine  current 
farm— A  Taiwanese  case  study 

M.Q.  Lee3,  C.N.  Lua,  H.S.  Huang  b* 

a  Department  of  Electrical  Engineering,  National  Sun  Yat-Sen  University,  No.  70,  Lienhai  Rd.,  Kaohsiung  804,  Taiwan,  ROC 
b  Department  of  Electrical  Engineering,  Ching  Yun  University,  No.  229,  Chien-Hsin  Rd.,Jung-Li  320,  Taiwan,  ROC 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  24  October  2008 
Accepted  21  January  2009 


Keywords: 

Reliability  analysis 
Cost  analysis 

Electricity  collection  systems 
Marine  current  farm 
Taiwan 


The  exploration  of  ocean  energy  for  electric  power  production  offers  a  sustainable  option  to  enhance  the 
use  of  renewable  energy.  In  this  article,  the  reliability  and  cost  analyses  results  of  electricity  collection 
systems  proposed  for  a  marine  current  farm  are  presented.  A  methodology  based  on  the  probability 
density  function  of  site  current  speed  is  developed  to  determine  the  speed  specifications  of  marine 
current  turbine.  Reliability  analyses  are  conducted  by  taking  electricity  collection  structure,  equipment 
failure  rate  and  probability  distribution  of  turbine  power  output  into  account.  Non-delivered  energy  cost 
in  conjunction  with  the  investment  cost,  power  loss,  operations  and  maintenance  costs  are  included  in 
the  cost  analyses.  Ocean  current  speed  data  measured  at  the  Taiwan  coastline  situated  in  the  Kuroshio 
stream  path  are  used  to  calculate  the  life-cycle  costs  of  the  studied  energy  collection  systems.  Simulation 
results  show  that  marine  turbine  parameters  can  be  effectively  specified,  and  a  sectionalized  radial 
collection  structure  provides  an  efficient  scheme  for  harnessing  ocean  energy. 
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1.  Introduction 

While  large  offshore  wind  farms  are  becoming  popular,  in  large 
areas  with  high  currents,  it  will  be  possible  to  install  marine 

*  Corresponding  author.  Tel.:  +886  3  4581196x5320;  fax:  +886  3  4594937.  turbines  in  groups  or  clusters  to  make  up  a  marine  current  farm. 

E-mail  address:  hshuang@cyu.edu.tw  (H.s.  Huang).  The  basic  physics  of  a  marine  current  turbine  (MCT)  are  similar  to 
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Nomenclature 

A  swept  area  of  the  marine  current  turbine  (m2) 

Q  initial  investment  cost  ($) 

Cp  power  coefficient 

CpF  present  value  cost  of  the  fixed  loss  cost  ($) 

(-pnde  present  value  cost  of  the  non-delivered  energy  cost 

(S) 

Cpom  present  value  cost  of  the  operation  and  main¬ 
tenance  cost ($) 

Cpv  present  value  cost  of  the  variable  loss  cost  ($) 

CF  capacity  factor 

a  critical  value  in  standard  table  calculated  based 
on  sample  size  n  and  significance  level  a 
f(v)  current  speed  probability  density  function 

FC  future  cost  being  paid  for  each  year  ($/year) 
irate  discount  rate 

LCC  life-cycle  cost  ($) 

rile  life-cycle  period  of  equipment  (year) 

Ns  number  of  operating  states 

NDE  non-delivered  energy  (MWh/year) 

Pe  electrical  power  output  (MW) 

Pe.avg  average  electrical  power  output  (MW) 

Pe<i  average  power  output  of  marine  current  turbine  i 
(MW) 

Pe<R  rated  electrical  power  output  (MW) 

PN  normalized  average  power 

P(x)  cumulative  normal  distribution  function 

Probs  probability  of  a  turbine  availability  state  (or  a 
generation  state) 

PW  present  worth  ($) 

Q.  forced  outage  rate  of  a  marine  current  turbine 

5„(x)  the  sample  cumulative  distribution  with  sample 
size  n 

TPI  turbine  performance  index 

Ui  unavailability  of  marine  current  turbine  i 

Vc  cut-in  speed  of  marine  current  turbine  (m/s) 

VF  cut-out  speed  of  marine  current  turbine  (m/s) 

VR  rated  speed  of  marine  current  turbine  (m/s) 

Greek  symbols 

standard  deviation  bounds 
rj  efficiency  of  drive  train 

/x  mean  value  of  the  normal  distribution  function  (m/ 

s) 

p  sea  water  density,  which  is  about  1024  (kg/m3) 

a  standard  deviation  of  the  normal  distribution 

function  (m/s) 


those  of  a  wind  turbine.  The  power  available  is  proportional  to  the 
working  fluid  density  and  cubic  of  its  speed.  MCTs  have  lower 
rotational  speeds  than  wind  turbines,  but  since  the  density  of  sea 
water  is  much  higher,  the  power  output  from  an  MCT  is  higher  than 
that  from  a  wind  turbine  of  equivalent  blade  size,  or  alternatively, 
MCTs  can  be  much  smaller  than  wind  turbines  to  generate  as  much 
power.  MCT  technologies  are  being  developed  in  many  countries 
[1].  Small-scale  systems  were  deployed  and  tested.  For  a  large- 
scale  system  which  covers  a  wider  area,  for  economic  reasons,  the 
current  speed  characteristics,  reliability  and  cost  of  electricity 
collection  system  have  to  be  investigated. 


The  power  output  of  an  MCT  depends  upon  the  mean  and 
distribution  of  current  speeds  at  the  site,  and  the  operation 
characteristics  of  the  MCT.  The  choice  of  cut-in,  cut-out  and  rated 
speeds  of  MCT  should  match  the  current  conditions  in  order  to 
bring  its  potential  into  full  play.  In  the  wind  farm  development,  it 
has  been  shown  that  wind  speed  at  a  location  can  be  statistically 
modeled  by  a  Weibull  probability  density  function.  In  [2], 
Jangamshetti  and  Rau  presented  a  method  of  matching  wind 
turbine  generators  to  site  wind  speed  characteristics.  This 
technique  identifies  optimum  turbine  speeds  from  a  turbine 
performance  index  (IP/)  curve  that  is  obtained  from  the  normal¬ 
ized  power  and  capacity  factor  curves. 

Establishing  a  marine  current  farm  requires  a  large  investment 
and  it  is  crucial  to  maximize  the  power  output  in  order  to  shorten 
the  pay  back  time.  Using  the  Kolmogorov-Smirnov  test  technique 
and  investigating  a  few  different  distributions,  it  is  concluded  that 
normal  distribution  is  more  adequate  for  describing  the  ocean 
current  speed  recorded  in  this  study.  This  difference  in  speed 
distribution  is  essential  and  should  be  paid  attention  if  techniques 
adopted  in  wind  turbine  design  are  to  be  used  in  the  MCT  design. 

Bahaj  and  Myers  proposed  a  method  to  assess  the  potential  of 
electricity  generation  at  Alderney  Race  in  the  Channel  Islands  [3]. 
This  work  quantified  the  resource  and  identified  currents  that 
could  be  used  for  the  establishment  of  various  size  arrays  in  the 
Race.  Predicted  resource  was  7.4  TWh  which  is  equivalent  to  2%  of 
UK  requirements  in  the  year  2000.  However,  in  that  work 
unavailability  and  power  loss  of  electricity  collection  scheme 
were  not  considered. 

Reliability  is  critical,  since  accessibility  to  the  offshore  marine 
current  farm  may  be  limited  during  certain  times  of  the  year  due  to 
adverse  weather  conditions.  Experiences  from  the  Middelgrunden 
offshore  wind  farm  reported  in  [4]  showed  that,  while  the 
availability  of  the  wind  turbines  was  around  98%,  the  availability 
of  the  whole  wind  farm  including  the  30-kV  collection  grid  was  at 
least  2%  lower. 

Sannino  et  al.  analyzed  the  availability  of  three  different 
collection  schemes  within  a  wind  farm  [5].  The  energy  not  produced 
due  to  outages  in  the  wind  farms  was  calculated  and  sensitivity  of 
the  results  to  variations  in  failure  rate  and  mean  time  to  repair  of 
cables  and  apparatus  was  presented.  In  [6],  a  comparative  analysis  of 
various  design  options  for  electricity  collection  system  of  large 
offshore  wind  farms  was  presented.  The  advantages  and  disadvan¬ 
tages  in  terms  of  their  steady-state  performance,  such  as  power 
losses,  voltage  level  changes,  contingency  analysis,  and  economics 
were  discussed.  Four  basic  designs  of  energy  collection  systems: 
radical  design,  single-sided  ring  design,  double-sided  ring  design 
and  star  design  were  analyzed.  According  to  the  steady-state  test 
results,  a  single-sided  ring  arrangement  is  a  better  option  for  a  large 
offshore  collector  system  as  it  achieves  fewer  losses  in  both,  normal 
and  contingency  operations.  However,  the  results  of  the  economic 
assessment  have  shown  that  the  cost  of  the  associated  redundant 
circuit  required  by  this  design  is  very  likely  to  become  an  economic 
barrier,  particularly  when  compared  with  the  cheaper  option  of  a 
typical  radical  design.  In  that  study,  the  economic  assessment  was 
based  only  on  the  investment  cost  and  system  power  loss. 
Unavailability  costs  were  not  considered. 

In  this  article,  the  reliability  and  cost  analyses  of  various 
electricity  collection  schemes  for  an  offshore  marine  current  farm 
are  presented.  In  the  analyses,  the  MCT  speed  parameters  are 
chosen  to  match  the  characteristics  of  offshore  current  at  northeast 
coast  of  Taiwan.  TPI  formulation  used  to  determine  suitable  turbine 
parameters  for  better  asset  utilization  is  derived.  With  the 
calculated  MCT  average  power  output,  life-cycle  cost  (LCC)  of  five 
different  electric  power  collection  schemes  are  presented.  The 
procedure  and  methodologies  involved  in  the  analyses  are  shown 
in  Fig.  1.  Details  are  given  in  the  following  sections. 
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Procedures  and  associated  methods  Data  outputs  served  as  input  of  the 

next  procedure  block 


Fig.  1.  The  procedure  used  for  reliability  and  cost  analyses  of  marine  current  farm  electricity  collection  systems. 


2.  Matching  turbine  specifications  with  current  characteristics 

2.2.  Identification  of  ocean  current  speed  distribution 

The  first  step  in  assessing  the  energy  potential  of  a  proposed 
marine  current  farm  is  to  identify  its  current  speed  distribution. 
The  hourly  ocean  current  speed  data  recorded  at  the  Su-ao  anchor 
station  are  used  in  this  study.  Fig.  2  shows  a  sample  distribution  of 
occurring  frequencies  of  the  current  speeds  at  30  m  below  sea 
level.  Higher  current  speed  locations  are  being  searched.  In  the 
following  presentation  the  data  shown  in  Fig.  2  are  used. 

Kolmogorov-Smirnov  test  (K-S  test),  which  quantifies  a 
distance  between  the  empirical  cumulative  distribution  function 
of  the  sample  and  the  cumulative  distribution  function  of  a 
hypothesized  distribution,  is  a  goodness-of-fit  test  and  appropriate 
for  testing  sample  data  against  a  continuous  distribution  [7].  In  this 
test,  if  the  value  of  the  maximum  deviation  exceeds  the  critical 
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Fig.  2.  Current  speed  distribution  of  Su-ao  station. 


value,  one  rejects  the  hypothesis  that  the  sample  is  from  the 
hypothesized  distribution.  It  is  used  in  this  study  and  distributions 
such  as  normal,  Weibull,  Rayleigh  and  other  distributions  were 
tested. 

If  the  hypothesized  distribution  is  a  normal  distribution,  the 
procedure  of  K-S  test  is  as  follows.  Given  a  sample  of  n  speed 
observations,  one  determines  D  =  maxx\P(x)  -  Sn(x)|,  where  S„(x)  is 
the  sample  cumulative  distribution  and  P(x)  is  the  cumulative 
normal  distribution  function  with  /x  =  X,  the  sample  mean,  and 
a2  =  s2,  the  sample  variance.  If  D  is  less  than  D“,  which  is  a  critical 
value  in  standard  table  calculated  based  on  sample  size  n  and 
significance  level  a ,  then  the  sample  data  follow  the  normal 
distribution.  If  D>D“,  one  rejects  the  hypothesis. 

For  the  1 009  data  shown  in  Fig.  2,  if  the  significance  level  is  0.05, 
then  =  D^0°059  =  0.0428.  D  =  maxx\P(x)  -  Sn(x)|  =  0.0252,  which 
is  less  than  D“.  Test  results  indicate  that  normal  distribution  is 
suitable  for  describing  the  characteristics  of  the  current  speeds 
measured  at  different  depth  of  the  studied  site. 

2.2.  MCT  speed  specification 

MCT  power  output  depends  upon  the  mean  current  speed,  cut- 
in,  cut-out  and  rated  speeds  of  MCT.  Thus,  there  is  a  need  for 
developing  a  method  to  identify  optimal  speed  parameters  of  MCTs 
in  order  to  deliver  higher  average  power  at  higher  capacity  factor. 
Table  1  shows  some  of  the  current  MCT  technologies  [8]. 

The  MCT  outputs  can  be  represented  by  their  V-P  (current 
velocity-power  output)  curves  as  shown  in  Fig.  3.  Fig.  3(a)  shows  a 
model  that  the  rated  speed  is  the  same  as  cut-out  speed.  Models 
MCT-A  to  MCT-G  shown  in  Table  1  have  such  V-P  curves.  MCT-H  is 
a  model  that  uses  rotor  pitch-control  and  has  a  flat  rated  output 
region  between  the  rated  speed  and  cut-out  speed  as  shown  in 
Fig.  3(b).  V-P  curve  shown  in  Fig.  3(b)  is  used  in  this  study  and  the 
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Table  1 

Parameters  of  some  marine  current  generators  [8]. 


Types 

Cut-in  speed 

(m/s) 

Rated  speed 
(m/s) 

Cut-out  speed 

(m/s) 

Rating 

(kW) 

A 

0.5 

2.58 

2.58 

6.95 

B 

1 

3.1 

3.1 

2000 

C 

0.7 

2.57 

2.57 

1520 

D 

0.7 

3 

3 

44 

E 

0.7 

2.3 

2.3 

1000 

F 

1.54 

3 

3 

400 

G 

0.7 

2.2 

2.2 

38 

H 

0.7 

3.1 

3.5 

2548 

power  output  can  be  expressed  as: 


On  the  other  hand,  if  the  CF  is  chosen  too  low,  then  although  the 
system  can  yield  a  higher  amount  of  output  at  rated  speed,  the 
turbine  would  seldom  operate  at  its  rated  capacity  and  lose  much 
of  the  energy  of  lower  speed  current. 

From  the  current  MCT  technologies  shown  in  Table  1 ,  one  could 
choose  an  MCT  with  the  highest  TPI  value  or  to  maximize  the  TPI  by 
specifying  the  speed  parameters  of  the  MCT.  In  this  study,  the  rated 
speed  vR  is  chosen  to  maximize  TPI  with  the  pre-specified  vc  and  vF, 
that  are  determined  based  on  site  current  speed  distribution  and 
the  Chebyschev’s  inequality  [9]  shown  in  the  following. 


fi\>Pcr)  < 


Pe(V) 


'0  forv<vcorv>vF 

1  , 

1  -pACprjV  fort'c  <v<vR 

1  ^ 

-pACpriVp  forvR<v  <vF 


where  p  is  the  sea  water  density,  A  is  the  swept  area  of  the  MCT,  Cp 
is  the  power  coefficient,  r\  is  the  efficiency  of  drive  train,  vc  is  the 
cut-in  speed  of  MCT,  vR  is  the  rated  speed  and  vF  is  the  cut-out 
speed. 

TPI  that  takes  both  capacity  factor  (CF)  and  power  output  into 
account  is  utilized  in  this  study  to  specify  MCT  speed  parameters.  It 
is  defined  as 


TPI 


PN  x  CF 


N.max 


X  CF 


max 


where  PN  is  the  normalized  average  power.  CF  indicates  average 
usage  of  the  assets,  is  equal  to  the  average  output  power  divided  by 
turbine’s  rated  output  power.  It  can  be  derived  based  on  the  speed 
distributions  determined.  If  cut-in,  cut-out  and  rated  speeds  are 
chosen  as  independent  variables,  and  f(v)  is  the  current  speed 
probability  density  function,  CF  is  expressed  as 


CF 


Pe.ave  _  Jo  Pe  (V)  f(V)  dv 
PeR  PeR 


J%(i/2)pACprii?f(v)dv  +  %(l/2)pACpTii$f(v)dv 

(\ /2)  pAC  pr]vl 


f(v)  dv 


where  pt  is  the  mean  value  of  the  normal  distribution  function  of 
current  speed,  a  is  the  standard  deviation  and  /3  is  the  standard 
deviation  bounds. 

Since  the  mean  value  and  standard  deviation  of  a  normal 
distribution  are  sufficient  to  make  statements  on  the  probability  of 
speeds  within  given  upper  and  lower  boundaries,  these  bounds  can 
be  used  to  specify  vc  and  vF. 

2.3.  Average  power  output  of  an  MCT  group 

The  wake  of  upstream  turbines  may  interfere  with  downstream 
turbines,  reducing  the  amount  of  energy  produced.  In  a  wind  farm, 
the  wake  effect  is  a  function  of  wind  turbines  spacing  (both 
downwind  and  crosswind),  wind  turbine  operating  characteristics, 
the  number  of  turbines,  size  of  the  wind  farm,  turbulence  intensity 
and  frequency  distribution  of  the  current  direction  [10].  In  a  marine 
current  farm,  distance  between  MCTs  may  be  large  in  order  to  avoid 
wake  effect.  To  simplify  the  reliability  model  of  a  large  marine 
current  farm,  the  wake  effect  is  neglected  and  a  group  of  MCTs  in  a 
neighboring  area  are  aggregated.  The  average  power  output  of  the 
group  is  calculated  by  using  the  capacity  outage  probability  table 
[  1 1  ].  As  the  name  suggests,  it  is  a  simple  array  of  capacity  levels  and 
the  associated  probabilities  of  existence.  For  a  group  of  m  MCTs,  each 
has  a  forced  outage  rate  of  Q,  the  probability  of  a  turbine  availability 
state  (or  a  generation  state),  Probs,  with  s  machines  up  and  m-s 
machines  down,  can  be  calculated  by  the  following  equation: 

s  m 

Probs  =  77(1  -  Q)  77  Q.  (6) 

1=1  j=s+l 


Details  of  the  CF  derivation  are  shown  in  Appendix  A.  A 
normalized  average  power  obtainable  from  MCT  is  defined  as  [2] 


P  e,ave 

( \/2)pACpri 


PeR  X  CF 
(1/2  )PACpT) 


Vp  x  CF 


The  average  power  output  of  the  MCT  group  is  the  summation 
of  multiplications  of  state  probability  and  amount  of  power  output 
in  that  state.  Repair  time  and  forced  outage  rate  could  be  high  due 
to  offshore  environment  and  affect  the  average  power  output. 


The  speed  specifications  of  marine  turbine  will  affect  CF  and 
they  should  be  determined  based  on  the  actual  current  speed 
characteristics.  If  a  higher  CF  is  desired,  then  turbine  rated  speed 
cannot  be  chosen  too  high  such  that  the  system  can  operate  at 
rated  power  output  condition  for  most  of  the  time,  but  in  this  case, 
much  of  energy  in  higher  speed  currents  would  not  be  harvested. 


3.  Reliability  and  cost  analyses  of  electricity  collection  schemes 

3.1.  Electricity  collection  schemes 

MCT  equipments  are  interconnected  in  strings  by  sub-sea 
cables.  The  strings  are  then  connected  to  the  platform  by  feeder 


Fig.  3.  V-P  curves  of  MCTs. 
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cables.  The  cables  together  with  necessary  switching  equipment, 
such  as  fuse,  breakers,  dis-connectors,  etc.,  form  the  collection  grid. 
There  are  various  possible  arrangements  for  the  energy  collection. 
Based  on  schemes  suggested  for  offshore  wind  farms  [5,6],  the 
authors  identified  five  basic  electricity  collection  designs  for  the 
studied  marine  current  farm. 

Radial  design  (scheme  a)  is  the  most  straightforward  arrange¬ 
ment.  Fig.  4(a)  shows  a  radial  design  with  circuit  breakers  at 
platform.  In  this  case,  the  MCT  equipment  shown  aggregates  five 
1  MW  MCTs.  The  maximum  number  of  MCT  equipments  on  each 
string  feeder  is  determined  by  the  capacity  of  the  generators  and 
the  maximum  rating  of  the  sub-sea  cable.  This  design  offers  the 
benefits  of  being  simple  and  inexpensive  because  the  total  cable 
length  is  smaller.  The  major  drawback  of  this  design  is  its  poor 
reliability  as  cable  or  switchgear  faults  at  the  platform  end  of  a 
radial  string  would  prevent  all  downstream  turbines  from 
exporting  power. 

Fig.  4(b)  illustrates  a  ring  design  (scheme  b),  where  two  strings 
are  interconnected  in  a  ring  system  with  a  few  sectionalizing 
breakers  and  the  ring  is  connected  to  the  platform  by  one  breaker. 
If  the  full  output  power  of  the  MCT  equipments  in  one  of  the  strings 
were  to  be  transferred  through  the  other  string  during  special 
situations,  then  the  cable  needs  to  be  sized  to  a  higher  rating. 

U-ring  design  (scheme  c)  shown  in  Fig.  4(c)  is  similar  to  the  ring 
design.  Each  string  of  MCT  equipments  is  connected  to  the 
platform  with  a  feeder  cable,  and  the  other  end  of  one  string  is 
connected  to  the  next  forming  a  ring  (U-ring).  The  redundancy 
offered  by  the  ringed  designs  is  beneficial  when  the  probability  of  a 


fault  and  the  associated  costs  are  higher  than  the  costs  associated 
with  additional  equipment,  especially,  when  the  repair  downtimes 
are  long.  However,  it  should  be  pointed  out  that  the  fault  current  of 
a  ring  structure  will  be  higher  than  that  of  a  radial  structure  and 
could  result  in  a  higher  requirement  in  switchgear  rating. 

The  star  design  (scheme  d)  illustrated  in  Fig.  4(d)  can  be  used  to 
reduce  cable  ratings  and  provide  a  high  level  of  security,  since  one 
cable  outage  only  affects  one  MCT  equipment.  However,  there  is 
additional  expense  in  longer  diagonal  cable  runs  and  some  short 
sections  of  higher  rated  cabling. 

The  fifth  scheme  will  be  presented  later.  In  Fig.  4,  a  three- 
winding  step-up  transformer  in  the  platform  is  also  shown.  This  is 
the  same  for  all  topologies  and  is  not  included  in  the  analysis 
presented. 

3.2.  Reliability  assessment  of  electricity  collection  system 

Reliability  analyses  can  be  conducted  with  relatively  simple 
calculations  by  using  the  network  reduction  method  [11],  which 
creates  a  sequence  of  equivalent  components  obtained  by  gradually 
combining  series  and  parallel  components.  Fault  tree  analysis  is  also 
a  popular  technique,  in  which  a  particular  failure  condition,  known  as 
the  top  event,  is  considered,  and  a  tree  is  constructed  by  identifying 
various  combinations  and  sequences  of  other  failures  that  lead  to  the 
failure  being  considered.  The  process  starts  by  first  identifying  a 
specific  mode  of  system  failure  to  be  assessed.  The  causes  of  this 
operation  mode  are  gradually  broken  down  into  an  increasing 
number  of  hierarchical  levels  until  a  level  is  reached  at  which  the 
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Fig.  4.  Electricity  collection  schemes. 
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behavior  or  effects  of  failure  of  the  basic  system  components  can  be 
identified.  The  fault  tree  analysis  is  adopted  in  this  work.  Fig.  5  shows 
a  fault  tree  that  MCT-A  cannot  deliver  power  to  the  platform  for 
scheme  a.  A  complete  model  of  each  collection  scheme  can  be  built 
based  on  the  interconnection  structure  of  the  scheme.  Combining  all 
the  basic  events  together  using  the  properties  of  Boolean  algebra  and 
the  logical  structure  of  the  fault  tree,  the  probability  of  the  top  event 
of  a  fault  tree  can  be  calculated. 

Different  reliability  indices  can  be  used  for  different  purposes, 
such  as  expected  generated  marine  energy,  average  service 
availability  index,  non-delivered  energy  (NDE),  and  average  service 
unavailability  index.  In  this  study,  the  NDE  index  defined  in  the 
following  is  adopted. 

NDE  =  X>e,A  (7) 

i= 1 

where  Pei  is  the  average  power  output  of  MCT  i,  D,-  is  its 
unavailability  and  Ns  is  the  number  of  operating  states. 

3.3.  Cost  analyses 

A  specific  design  choice  may  have  a  complex  effect  on  the  project 
financial  performance,  affecting  capital  costs,  taxes,  insurance, 


energy  revenue,  maintenance  costs,  and  government  subsidies.  In 
this  study  a  simplified  calculation  based  on  minimizing  LCC  of 
ownership  and  operation  of  the  network,  is  used  to  compare 
different  design  options.  A  scheme  with  minimum  LCC  that  balances 
the  capital  investment  against  the  costs  of  system  operation  and 
maintenance  (O&M)  and  energy  losses  is  preferred. 

There  are  three  unique  types  of  energy  losses,  each  affecting  the 
revenue  model  in  a  different  way,  they  are  [12]: 

(1)  Fixed  losses  which  do  not  vary  with  turbine  output.  These 
losses  are  primarily  transformer  excitation  losses. 

(2)  Variable  load  losses  which  are  ohmic  losses  in  cables  and 
transformers. 

(3)  Energy  non-delivered  due  to  a  constraint  imposed  by  electrical 
system  unavailability. 

Present  values  of  the  costs  and  losses  in  the  system  life  time 
translate  annual  costs  and  lost  energy  production  into  an 
equivalent  initial  capital  cost.  The  LCC  of  a  studied  system  can 
be  defined  as 

LCC  =  Cj  +  CpoM  +  CpF  +  CpV  +  CpjvDE  (8) 

where  C7  is  the  initial  investment  cost,  CPOm  is  present  value  cost  of 
the  operation  and  maintenance  cost,  CPF  is  present  value  cost  of  the 
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Fig.  5.  A  fault  tree  example  for  scheme  a. 


fixed  loss  cost,  CPV  is  present  value  cost  of  the  variable  loss  cost  and 
Cpnde  is  present  value  cost  of  the  NDE  cost. 

The  cost  transformation  from  future  value  to  present  value  is 
based  on  following  equation 


PW 


PC 


(1  +  irate)nk  -  1 
irate  x  (1  +  irate)nlc 


where  PW  is  the  present  worth,  irate  is  the  discount  rate,  nlc  is  the 
life-cycle  period  of  equipment  and  PC  is  the  future  cost  being  paid 
for  each  year. 


4.  Analyses  results  and  discussions 


4.1.  MCT  speed  specification  results 

Ocean  current  data  shown  in  Fig.  2  are  used  for  turbine  speed 
parameter  specification.  As  shown  in  Fig.  6,  if  the  probability  of  the 
speed  of  ocean  current  exploited  to  generate  electric  power  is 
chosen  to  be  0.8,  then  the  probabilities  outside  the  upper  and  lower 


bounds  are  0.1  and  0.1,  respectively,  and  from  Chebyschev’s 
inequality  shown  in  (5),  l/(2/32)  =  0.1,  /3  is  2.236.  The  speeds 
corresponding  to  the  lower  and  upper  bounds  are  pt  -  /3cr  =  0.301 5 
(m/s)  and  pt  +  /3cr=  1.4910  (m/s)  respectively.  These  could  be 
chosen  as  the  cut-in  and  cut-out  speeds  of  the  turbine.  However, 
after  considering  the  current  marine  turbine  technology  (Table  1), 
the  cut-in  speed  is  set  to  0.5  m/s. 

With  the  specified  vc  and  vF ,  the  CP,  PN  and  TPI  curves  are 
calculated  and  shown  in  Fig.  7.  It  can  be  seen  that  PN  and  CF  curves 
reach  their  peaks  at  different  rated  speeds  and  implying  that  a 
compromise  should  be  made  to  achieve  the  best  result.  Since  TPI 
curve  takes  CF  and  normalized  power  factors  into  to  account,  its 
peak  can  be  used  to  specify  the  rated  speed  of  the  turbine.  The 
turbine  specifications  are  summarized  in  Table  2. 

4.2.  Average  power  output  of  a  MCT  group 

The  MCT  group  shown  in  Fig.  8  has  five  1  MW  MCTs  and  each 
has  a  forced  outage  rate  of  0.05,  excluding  the  consideration  of 
bus  bar  and  transformer  failures  for  simplicity,  the  probabilities 


Fig.  6.  Normal  distribution  representation  of  current  speeds. 


Fig.  7.  The  capacitor  factor,  normalized  power,  and  TPI  curves. 
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Table  2 

MCT  specifications. 


Parameters  and  calculated  results 

Values 

Cut-in  speed  (m/s) 

0.5000 

Rated  speed  (m/s) 

1.0100 

Cut-out  speed  (m/s) 

1.4910 

Normalized  power 

0.6548 

Capacity  factor  (CF) 

0.6355 

Prated=1  MW 


MCT 

Equipment 


of  output  states  and  average  power  output  of  the  group  are 
shown  in  Table  3.  The  calculated  CF  shown  in  Table  2  is  used 
to  calculate  the  turbine  average  power  output.  If  the  outage 
rate  of  each  MCT  is  neglected,  then  the  average  power  output 
of  the  MCT  group  would  be  0.1589  MW  higher  which  results 
in  an  overestimate  of  energy  yield  by  1391.964  MWh  per 
year. 

4.3.  Reliability  analyses  results 

The  reliability  indices  of  the  four  collection  schemes  described 
in  Section  3.1  are  calculated  using  the  failure  rate  and  repair  time 
assumptions  shown  in  Table  4. 

Using  the  fault  tree  technique,  the  NDEs  of  the  four  collection 
schemes  are  calculated  for  cases  with  and  without  fuses  which 
are  installed  between  the  MCT  equipment  and  feeder.  It  can  be 
seen  from  Table  5  that  fuses  with  a  small  additional  investment 
can  be  used  to  enhance  the  reliability  of  the  energy  collection 
system.  In  the  following  discussion,  it  is  assumed  that  MCT 
equipment  is  always  equipped  with  a  fuse.  As  can  be  seen 
scheme  c  (U-ring  design)  has  the  best  reliability  performance. 
Compare  to  scheme  a  (radical  design),  its  NDE  is  2749  MWh/year 
lower.  The  results  indicate  that  when  the  repair  time  is  long, 
NDE  would  provide  useful  information  in  choosing  an  effective 
electricity  collection  scheme.  However,  it  should  be  pointed  out 
that  the  fault  currents  of  ring  structures  will  be  higher  than  that 
of  radial  structure  and  could  result  in  higher  requirement  in 
switchgear  rating. 


Table  4 

Reliability  data  of  the  components  involved. 


Components 

Failure  rate  (fr/yr) 

Repair  time  (h) 

Feeder 

0.015 

2160 

FeederMCr-MCT 

0.01 

2160 

CableB, c„ .  ,h 

0.01 

2160 

MV  breaker 

0.025 

240 

Platform  breaker 

0.025 

72 

MCT  equipment 

0.1 

168 

Table  5 

NDEs  of  different  electricity  collection  schemes. 

Scheme 

NDE  (MWh/year) 

NDE  (MWh/year) 

(without  fuse) 

(with  fuse) 

(a) 

7437 

4701 

(b) 

3824 

3433 

(c) 

2343 

1952 

(d) 

3447 

3104 

4.4.  Cost  analyses  results 

Connection  structures  with  various  cable  sizes  shown  in  Fig.  4 
are  studied.  The  lengths  of  submarine  cables  are  1.5  km  between 
two  adjacent  MCT  equipments  and  5  km  between  MCF  and 
platform  (the  800  mm2  cable).  Information  of  XLPE  cable 
impedances  and  shunt  charging  are  available  in  [13].  The  costs 
of  equipment  involved  are  listed  in  Table  6  [14].  Cost  analyses  are 
conducted  based  on  the  following  assumptions:  electric  price  is 
0.15  $/l<Wh,  O&M  cost  is  0.3  M$/year,  discount  rate  is  5%  and  the 
economic  life  is  20  years.  If  the  fixed  loss  is  neglected,  Table  7 
shows  the  costs  and  LCC  calculated.  Variable  loss  cost  is  computed 
based  on  line  losses  obtained  from  power  flow  solutions  of  various 
electricity  grid  schemes.  From  Table  7,  it  can  be  seen  that  scheme  c 
(U-ring  design)  has  the  lowest  LCC.  The  U-ring  design  has  the 
highest  investment  cost,  but  lowest  variable  loss  and  NDE  costs.  On 
the  other  hand,  radical  design  (scheme  a)  has  the  lowest 
investment  cost,  but  higher  variable  loss  and  NDE  costs.  Since 
NDE  costs  are  higher  than  variable  loss  costs  and  contribute  11- 


Table  6 

Costs  for  equipment  involved  [14]. 


Equipment 

Cost  ($) 

Breaker 

135000 

Fuse 

10000 

95  mm2  conductor 

152/m 

150  mm2  conductor 

228/m 

240  mm2  conductor 

250/m 

300  mm2  conductor 

300/m 

400  mm2  conductor 

381/m 

630  mm2  conductor 

571/m 

800  mm2  conductor 

600/m 

Table  3 

Probability  of  the  turbine  availability  and  average  power  output  of  the  MCT  group  shown  in  Fig.  8. 


State 

No.  of  turbines 
in  service 

Total  rated 
capacity  (C,  MW) 

State  probability  Ps 

Average  power  output 
(MW)  (C  x  Ps  x  CF) 

0 

0 

0 

3.125E— 07 

0 

1 

1 

1 

2.96875E— 05 

1.88664E— 05 

2 

2 

2 

0.001128125 

0.001433847 

3 

3 

3 

0.021434375 

0.040864636 

4 

4 

4 

0.203626563 

0.517618723 

5 

5 

5 

0.773780938 

2.45868893 

Total 

1.000000000 

3.018625003 

2020 
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Table  7 

Costs  of  different  collection  schemes. 


Scheme 

Investment  costs  (M$) 

Variable  loss  costs 

Not-delivered 

costs 

energy 

O&M  costs 

Life-cycle  costs  LCC  (M$) 

(M$/year) 

PW  (M$) 

(M$/year) 

PW  (M$) 

(M$/year) 

PW  (M$) 

(a) 

17.99 

0.71 

8.848 

0.71 

8.848 

0.3 

3.739 

39.425 

(b) 

19.78 

0.33 

4.113 

0.52 

6.480 

0.3 

3.739 

34.112 

(c) 

24.87 

0.068 

0.847 

0.29 

3.614 

0.3 

3.739 

33.070 

(d) 

22.45 

0.200 

2.492 

0.47 

5.857 

0.3 

3.739 

34.538 

5  MW 


ABC 


D  MCT 


Fig.  9.  Collection  scheme  (p). 


22%  (present  worth)  to  LCC,  it  will  affect  the  decision  in  adopting  an 
electricity  collection  system  for  marine  energy  farm  development. 

Another  collection  scheme  (scheme  p)  modified  from  scheme  c 
by  opening  the  loop  and  resizing  the  conductor  is  shown  in  Fig.  9. 
The  LCC  of  this  scheme  is  31.727  M$,  which  is  lower  than  those  of 
other  schemes.  Scheme  p  has  a  better  reliability  performance  when 
compared  with  that  of  scheme  a.  This  scheme  also  has  a  lower 
investment  cost  than  that  of  scheme  c.  Therefore,  scheme  p  is 
capable  of  keeping  a  good  balance  between  cable  cost  and  system 
reliability.  Using  the  techniques  presented,  one  could  pursue 
better  design  with  lower  LCC  for  marine  energy  collection. 


The  capacity  factor,  CF 

_Pe,ave  _fo  Pe{V)f(V)  dv 
PeR  PeR 

C  ( V2)/qACpi? V3  m  dv  +/;;(!  /2)pACptlvl  f(v)  dv 

(1/2)  pACprjv3 


1  CVr 

/  v3f(v)dv  + 
V3RJvc 


f(v)  dv  = 

(v  ~  ^)2 

2a2 


CF!  +  CF 
^  dv 


2 


(A.2) 

(A.3) 


5.  Conclusion 


Let  x  =  (v  -  ii)/V2cr,  and  then  dx  =  dv/V2cr 


A  brief  review  of  offshore  wind  farm  design  is  given  and  an 
analytical  procedure  for  cost  analyses  of  electricity  collection 
systems  with  emphasis  on  system  reliability  of  a  marine  current 
farm  is  presented.  A  new  methodology  is  proposed  to  determine  the 
speed  specifications  of  MCT.  Analytical  results  show  that  the 
probability  density  function  of  ocean  current  speeds  used  in  this 
study  is  a  normal  distribution,  and  based  on  this  distribution  TPI 
formulation  is  derived  and  used  to  determine  turbine  speed 
specifications  in  order  to  achieve  higher  energy  output  and  capacity 
factor.  The  importance  of  unavailability  cost  is  highlighted  in  the  test 
result.  It  shows  that,  among  the  collection  schemes  investigated,  U- 
ring  design  has  the  best  reliability  performance  and  lower  LCC. 
However,  numerical  result  also  shows  that  a  simple  modified  scheme 
could  also  provide  an  effective  solution  for  marine  energy  collection. 

Appendix  A 


CF 
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where  a  =  (vc  -  fi)/V2a,  b  =  (vR  -  p,)/V2a 
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According  to  Steen’s  deduction  for  integration  Jq  exp(-x2) 
/(x)dx  by  using  Gaussian  quadratures  in  1969  [15],  CF  can  be 
calculated  by  the  following  equation: 


Normal  distribution  is  expressed  as: 
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where 

/0(a)  = 
/1(a)  = 
/2(a)  = 
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erf  (a)  is  the  error  function,  it  is  defined  as  the  following 
equation: 


2  /*“ 

er/(a)=— =/  exp(-r)dt 
do 


(A.l  1 ) 


In  a  similar  way: 

Let  x  =  (v  -  /a)/V2<j,  and  then  dx  =  dv/V2cr 


CF7  = 


f  e(-x2)dx 

2tt<j  Jb 


V2j to 

where  h  is  defined  as  before,  c  =  (vF  -  /jl)/V2o 

CF2  =  f=(I0(c)  -  10(b)) 


(A.12) 


(A.13) 


(A.14) 


From  the  above  deduction,  it  indicates  that  CF  is  a  function  of  a, 
b  and  c.  Therefore,  CF  is  also  a  function  of  cut-in  speed  vc ,  rated 


speed  Vk  and  cut-out  speed  vF,  and  can  be  described  by  a  closed 
form  function. 
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